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ABSTRACT 

We present high time-resolution multicolour observations of the quiescent soft X- 
ray transient V404 Cyg obtained with ULTRACAM . Superimposed on the secondary 
star's ellipsoidal modulation are large flares on timescales of a few hours, as well as 
several distinct rapid flares on timescales of tens of mins. The rapid flares, most of 
which show further variability and unresolved peaks, cover shorter timescales than 
those reported in previous observations. The power density spectrum (PDS) of the 
5 s time-resolution data shows a quasi-periodic oscillation (QPO) feature at 0.78 mHz 
(=21.5 min). Assuming this periodicity represents the Keplerian period at the transi- 
tion between the thin and advective disc regions, we determine the transition radius. 
We discuss the possible origins for the QPO feature in the context of the advection- 
dominated accretion flow model. 

We determine the colour of the large flares and find that the i'band flux per 
unit frequency interval is larger than that in the g' band. The colour is consistent 
with optically thin gas with a temperature of ~8000K arising from a region with an 
equivalent blackbody radius of at least 2 Rq , which covers 3 percent of the accretion 
disc's surface. Our timing and spectral analysis results support the idea that the rapid 
flares (i.e. the QPO feature) most likely arise from regions near the transition radius. 
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1 INTRODUCTION 

Soft X-ray transients (SXTs) are a subclass of low-mass X- 
ray binaries (LMXBs) that are characterized by episodic, 
dramatic X-ray and optical outbursts, which usually last for 
several months and recur on a timescale of decades. In the in- 
terim the SXTs are in a state of quiescence during which the 
optical emission is dominated by the luminosity of the faint 
companion star jvan Paradiis fc McClintock^ll995^ . In qui- 
escence the optical lightcurves exhibit the classical double- 
humped ellipsoidal modulation, which is due to the differing 
aspects that the tidally distorted secondary star presents to 
the observer throughout its orbit (Shahbaz ct al. 2003). 

In outburst the high/soft and the low/hard X-ray states 
are commonly seen. In the high/soft state, X-ray emission 
is dominated by thermal emission from an accretion disc ex- 
tending close to the last stable orbit around a black hole. In 
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the low/hard state, the inner disc is believed to be truncated 
and emission appears to arise from an extended corona. Sim- 
ilar ideas are involve d for the advective models proposed fo r 
the quiescent state iNaravan. Garcia fc McClintockll200ir) . 
but with the disc truncated at larger radii. The states of 
SXTs are also classified by their X-ray timing properties; 
the high/soft state shows low level red noise, whereas the 
low/hard state and very high states exhibit band-limited 
noise but with a low-frequency bre ak at ~0. 02-30 Hz^ and 
sometimes superposed QPOs (see Ivan der KlisI Il995l and 
tWiinands & van der Klis 19^. 

Although the spectral and timing properties during the 
quiescent state have not been well studied, one might expect 
similar properties, as the structure of the flow is believed 
to be similar to that in the low/hard state. The relatively 
high quiescent X-ray luminosity of V404 Cyg with respect 
to other SXTs indicates that there is still some continuing 
accretion. This is also supported by the short-term variabil- 
ity seen at all wavelengths in several quiescent photometric 
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Table 1. Log of ULTRACAM observations. 



Date 


Run # 


# Pts 


Exposure 


JD 2452520.5 + 


Orbital phase 


n 
u 


i 


Comments 








Time (s) 


start - 


end 


mid 








nn /n /no 

uy/y/uz 


run 09(a) 


1200 


5.0 


26.96793 - 


27.05119 


0.855 


19.836 


16.570 


Poor seeing at end of night 


10/9/02 


run 10(a) 


1266 


5.0 


27.87320 - 


27.94646 


0.993 


19.755 


16.615 




10/9/02 


run 10(b) 


1152 


5.0 


27.97038 - 


28.03704 


0.159 


19.920 


16.658 


Good seeing 


12/9/02 


run 12(a) 


1192 


5.0 


29.85483 - 


29.92379 


0.293 


19.545 


16.406 




12/9/02 


run 12(b) 


1282 


5.0 


29.92406 - 


29.99824 


0.310 


19.641 


16.467 




12/9/02 


run 12(c) 


848 


5.0 


29.99853 - 


30.04770 


0.316 


19.695 


16.516 




13/9/02 


run 13(a) 


1082 


5.0 


30.83847 - 


30.90107 


0.451 


19.968 


16.690 




13/9/02 


run 13(b) 


948 


5.0 


30.90135 - 


30.92339 


0.460 


19.879 


16.656 




14/9/02 


run 14(a) 


8544 


0.237 


31.83812 - 


31.86152 


0.602 




16.768 




14/9/02 


run 14(b) 


2791 


0.237 


31.86162 - 


31.86902 


0.605 




16.669 


Noisy data, not used in PDS 


14/9/02 


run 14(c) 


9155 


0.237 


31.86936 - 


31.89443 


0.607 




16.649 




14/9/02 


run 14(d) 


7236 


0.237 


31.89942 - 


31.91924 


0.611 




16.754 




14/9/02 


run 14(c) 


8773 


0.237 


31.91932 - 


31.94334 


0.615 




16.708 





studies of V404Cyg dPavlenko et al.lll996l : Zurita, Casares 
& Shahbaz 2003). Sub-orbital variability with an amplitude 
of 0.10-0.20 mag 1, Wagner et al. ,1994) which is part l y due 
to a ~6hr quasi-per iodic oscillation dCasares et alJ Il993l: 
iPaylenko et alj Il996l) is present in the optical. The Ha 
line profile also changes significantly on short timescaies 
iCasares fc Charles 1992) and seems to be correlated with 
the contirmum (Hynes et al. 2002.). Short timescale v aria- 
tions also appear in the infrared dSanwal et alj[l'996t) . At 
X-ray energies, ROSAT (0. 1-2.4 keV) saw changes of a fac- 
tor of 10 on timescaies <0.5 day. More recent observations 
with Chandra (0.3-7 keV) s howed variability of a factor of 2 
in a few ksec jKong et aljEb02l) . Finally V404Cyg is also 
a variable radio source on times caies of days, wi th typical 
variability in flux of 0.1-0.8 mjy jHiellminej|2000h . 



The origin of the variability remains uncertain, the most 
likely explanations are magnetic reconnection events in the 
disc or optical emission from an advective region. The other 
possibilities such as reprocessed X-ray variability and flick- 
ering from the accr etion stream impact point are less likely 
JZurita et al.ll2003ll . To determine if just one type of vari- 
ability dominates or if the variability is a combination of 
mechanisms on different timescaies, it is therefore impor- 
tant to perform a comparative study of the class. 



In the optical, the faintness of the SXTs has usually 
limited the time resolution one can practically achieve. Fur- 
thermore, observations have normally been constrained to 
a single waveband. ULTRACAM with 4-m class telescopes 
provides a unique opportunity not only to study photometric 
variability on the shortest timescaies, but also in more than 
one waveband simultaneously. Therefore we have embarked 
upon a programme to study the fast multicolour variability 
of quiescent SXTs with ULTRACAM . Our main aims are to 
study the multicolour variability of the brighter SXTs and to 
study variability in the fainter ones, faster than previously 
possible. We present here high time resolution multicolour 
observations of one of the brightest SXTs, V404 Cyg. 



2 OBSERVATIONS AND DATA REDUCTION 

Multicolour photometric observations of V404Cyg were 
taken with ULTRACAM on the William Herschel Telescope 
atop La Palma during the period 2002 Sep 9 to 13. UL- 
TRACAM is an ultra-fast, triple-beam CCD camera. The 
light is split into three wavelengths colours (blue, green and 
red) by two dichroic beamsplitters and then passes through 
a filter. The detectors are three back-illuminated, thinned, 
Marconi frame-transfer 1024x1024 active area CCD chips 
with a pixel scale of 0.3 arcsecs/pixel. The CCDs are cooled 
using a three-stage peltier device and water chiller, resulting 
in negligible dark current (especially at these short integra- 
tion times). With the frame transfer mo de, the dead-time 
is ess entially zero (for further details see iDhillon fc MarshI 
l200lD . 

Our observations were taken using the Sloan u' , g' and 
i' filters with effective wavelengths of 3550 A , 4750 A and 
7650 A respectively. For the first 5 nights (Sep 9 to 12), we 
used an exposure time of 5.0 sec, which was short enough to 
give reasonable counts in the g' and i' bands. Given the faint- 
ness of the object (_B=20.6, y=18.4) and the short exposure 
times, few counts were obtained in the u' band. On the last 
night (Sep 13) we decreased the exposure time to 0.237 sec. 
Only data in the i' band were usable. From hereafter, we will 
refer to the data taken with exposure times of 5.0 sec and 
0.237 sec as the "slow" and "fast" data respectively. 

The weather conditions during the observing run were 
on the whole very good. The night of Sep 9 sufTered from 
severe seeing (2") during the end of the night. No useful 
data could be obtained on Sep 11 due to cloud. The night of 
Sep 10 was p hotometric, so we observed the fiux standard 
BD-l-35 3659 llSmith et alJl2002l) in order to determine the 
u , and i' band instrumental zeropoints, which were then 
used to calibrate the data. On Sep 14, due to an earthing 
problem, the electronic bias level occasionally increased by 
a factor of 2, however, this only affected runs 14(b) and (c). 
A log of the observations is given in Table0 

The pipeline reduction procedures were used to debias 
and fiat-field the data. The same pipeline was also used 
to determine lightcurves for V404 Cyg and several compari- 
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son stars by extracting the counts using aperture photome- 
try. As there is a hne-of-sight contaminating star l.Sarcsec 
North of V404 Cyg, we obtained the combined counts using 
an aperture sufficient enough to encompass both stars. We 
determined hghtcurves with different sized apertures, rang- 
ing from to f.8 to 4.2arcsec. The optimal and most rehable 
results, given the contaminating star, were obtained using a 
large aperture of 3 arcsec. The count ratio of V404 Cyg with 
respect to the local standard (3.6 " S 75.8 " E of V404 Cyg) 
was then determined by subtracting the count ratio of the 
contaminating star with respect to the local standard de- 
termined from images taken under good seeing conditions 
(0.7arsecs), from the combined count ratio of V404Cyg 
(i.e. V404 Cyg -I- line-of-sight star) with respect to the lo- 
cal standard. The local standard was chosen to have a sim- 
ilar colour to the contaminating star. The count ratio of 
the contaminating star to the local standard is 0.065 and 
0.143 in the and i' bands respectively. The magnitude of 
V404 Cyg was then obtained using the observed magnitude 
of the local standard, g' =16.90 and i' =15.12 (<0.01 mag 
uncertainty). The magnitude of the contaminating star was 
(?' =19.87 and i' =17.23. As a check of the photometry and 
systematics in the reduction and extraction procedures, we 
also computed lightcurves of a comparison star (13.2 "S, 
17.1 "E of V404Cyg). The magnitude of the comparison 
star was g' =18.90±0.03, i' =16.55±0.02. The mean g' and 
i' band magnitudes of V404 Cyg were 19.77 and 16.57 respec- 
tively. We estimate the relative photometric accuracy to be 
3.0 and 0.9 percent for the 'slow' g' and i' band respectively 
and 6.0 percent for the "fast" i' band data. 

By co-adding all the images from run 10(b) on 
Sep 10 taken under good seeing conditions, we estimate 
?t'=22.6±0.7 for V404Cyg. 



3 BACKGROUND TO V404 Cyg 

The X-ray transient GS 2023-1-338 was discoverd during its 
ourburst in May 1989 by the All Sky Monitor aboard the 
Ginga satelite (.Making .198911 . Its high X-ray luminonsity 
suggested that, like other X-ray transients, the system was 
a close binary with an accreting compact object. The X- 
ray source was soon identifed with the known variable star 
V404 Cyg , which had bee n classified as a nova after its 1938 
outburst iMarsdenlll989l) . 

Spectroscopic observations taken when V404 Cyg was 
in quiescence revealed the companion to the X-ray source 
to be a G or early K-type star in an orbit with a pe- 
riod of 6.473 d, which when combined with the secondary 
star's radia l velocity curve implies a ma ss function of 
6.3±0.3Mr7. llCasares. Charles fc Navlo^ll992^ . Later on, ad- 
ditional measurements determined the secondary star's ro- 
tational broadening, and refined the spectral type of the sec- 
ondary star and the orbital parameters; the orbital period 
is Porb=6.4714 d and semi-amplitude of the secondary star's 
radial velocity curve is 7^2= 208.5 km s~^, giving amas s func- 
tion of f(M)=6.08±0.08 M,7. llCasares fc Charleijll994^ . Since 
this mass function, which provides a firm lower limit to the 
compact object's ma ss, exceeds the maximum mass of a neu- 
tron star f^3.3M^: rCook. Shapiro fc Teukolskv^ll994^ ■ the 
compact object in V404 Cyg is almost certainly a black hole. 

Optical photometry of V404 Cyg in the /-band revealed 
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Figure 1. The observed g' and i' flux lightcurves of V404 Cyg 
and the comparison star. The solid line is the secondary star's 
ellipsoidal modulation (see text). 
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the secondary star's ellipsoidal variation iWagner et alJ 
Il992h but these were, however, severely contaminated by 
short-ter m variability. Therefore , in an effort to avoid these 
problems, IShahbaz et alJ lll994ll observed the lightcurve of 
V404 Cyg in the isT-band , where the flux from the KOIV star 
is substantially greater iShahbaz et aljri996r) . The analysis 
of the ellipisoidal variations, which results from the tidal 
disortion of the secondary star gives the binary inclination 
angle, combined with spectroscopic observations yields the 
mass of the binary components. They concluded that the 
contamination in the infrared was small and deduced the 
blac k hole mass to be 12 ±2 Mq. 

ICasares et alj (^g^) estimated the interstellar redden- 
ing to V404Cyg, by comparing its observed colours with 
those of a KOIII star, after allowing for the accretion disc 
contribution to the observed flux. They obtained Ay ~ 4 
and a lower limit of Ay ~ 3.6. However, as one can see, 
the reddening determined in this way depends heavily on 
the disc contamination in the optical, which the authors 
measure with large uncertainites. A more accurate deter- 
minat ion of the reddening was obtained by IShahbaz et alJ 
il994) who obtained the JT-band lightcurve of V404 Cyg in 
quiescence and determined the binary inclination angle by 
modeling the secondary star's ellipsoidal modulatio n. As- 
sumi ng that the secondary star is a stripped giant iKind 
Il993h they also determined limits to the reddening and dis- 
tance to V404 Cyg, by matching the model stripped giant 
luminosity to that of the secondary star as predicted by the 
ellipsoidal model. Assuming no disc contamination in the K- 
band, a reasonable assumption given that the disc contam- 
ination measured in the iC'-band was found to be negligible 
iShahbaz et aLlFiggdl . and 10 percent in the l/-band, they 
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Figure 2. In the top panel we show the flare flux density Ughtcurve in the g' (dashed line) and i' (solid line) bands, obtained by 
subtracting the reddened secondary star's ellipsoidal modulation from the observed lightcurves. For clarity the flare lightcurves have 
been re-binned to a time resolution 1 min. The uncertainties in the g' and i' lightcurves are 0.02 mjy and 0.01 mjy respectively. The 
letters mark the rapid flare events used to determine the flare properties in Tablel^ The middle panels show the flux density ratio fgi / f^i 
and the bottom panel shows the projected blackbody radius of the region producing the flares. 



obtained 2.2< Ay <3.3 (90 percent confidence). Although 
the determination of tiie reddening tiirougii this method is 
model dependant, tlie uncertainites in the disc fraction and 
model dependant parameters such as the binary masses and 
inclination angle are less in the infrared compared to in the 
optical. Given the uncertainties in determining the redden- 
ing using optical and infrared data, throughout this paper 
we will use a range of values for the reddening, 2.2< Ay <4.0 
and a central value of Ay =2.8. 



4 THE LIGHTCURVE OF V404 Cyg 

The optical lightcurve of V404 Cyg is dominated by the sec- 
ondary star's ellipsoidal modulation, which is due to the 
observer seeing differing aspects of the tidally distorted sec- 
ondary star. This modulation gives rise to two maxima and 
minima per orbital cycle. However, superimposed on the el- 
lipsoidal modulation are the large 6 hr flares seen previously 



jPavlenko et aljll99d) and also many rapid flares (~0.5hr), 
similar to what is observed in the short orbital period SXTs 
JZurita et al.ll2003l: iHvnes et aLl liboS). It is clear that the 
secondary star's ellipsoidal modulation is not linked to the 
variability, therefore if we want to determine the flux of 
the flares, its contribution must be first removed from the 
lightcurves. 

Figure 1 shows the observed data. Large amplitude 
flares which last a few hours or more, are observed only 
on Sep 10, 12 and 13 and are most likely relate d to the well 
known 6hr flare/QPO ijC asarcs et al. 1993; Pa vlenko et alJ 
Il996l: iHvnes et aLn2002ir On Sep 9 no large flare is seen. 
In order to isolate the flux of the flare we assume that the 
light produced by a flare is simply added to the quiescent 
spectrum, where the quiescent spectrum is the light from the 
secondary star's ellipsoidal modulation and the non-variable 
accretion disc. 

We first de-redden the observed magnitudes using the 
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central value for Ay of 2.8 (se e section|^ and the ratio 
Ay/E{B - 1/)=3.1 JSeatorll979l '). giving g' and i' extinction 
values of 3.34 and 1.79 mags respectively and then convert 
the Sloan AB magnitudes to flux density . We u se the X-ray 
binary model described in lShahbaz et al.1 ll2003|) with the bi- 
nary m ass ratio and inclination a ngle given in |Casares et alJ 
il993h and lShahbaz et al] ||l99j) respectively, to determine 
the ellipsoidal modulation. The model is not set up to com- ^ 
pute Sloan magnitudes, so we first compute the Johnson- 
Cousins magnitudes and then co nvert to Sloan magnit udes B 
using the colour transformations iFukugita et alj[l996ll . We 
shift the model ellipsoidal lightcurve so as to fit the lower- ^ 
envelope of the observed lightcurve. We use the data at 
phase 0.85 (Sep 9) to anchor the fit, since on this night 
the large 6hr flare is not seen and so the lower-envelope 
of these data most likely represents the light from the sec- 
ondary star and non- variable accretion disc. Figure 1 shows 
the data and the scaled ellipsoidal model. Clearly the model 
fails to fit the g' and i' data near phase 0.0 (Sep 10). This 
suggests that the accretion disc light is variable, most likely 
due to a residual superhump modulation as is observed in 
other quiescent S XTs (sec Lcibowitz, Hcmar & Orio 19^; 
IZurita et al.ll20o3) . This will obviously introduce uncertain- 
ties in the subtraction procedure, however, we can obtain a 
good estimate for the flux of the flares, since on Sep 12 and 
13 the model fit partially matches the end and beginning 
of the large fiares respectively. Finally, the scaled ellipsoidal 
modulation is subtracted from the de-reddened lightcurve 
yielding the colour of the fiares (see Figure 2). 

One should note that the conversion from magnitudes 
to fiux density is colour dependant. In order to estimate 
the scale of this effect, we determine the fiuxes obtained 
using spectra with different power-law indices (F^, oc u"). 
The spectrum of the fiares has a mean flux density ratio of 
fg' / fi' ^^0.40 (see Figure 2) which corresponds to a spec- 
trum with a power-law index of ~-2. The ellipsoidal model 
has a colour (<?'-«' ) = 1.67 which corresponds to a spectrum 
with a power-law index of -2.8. Using these models, we esti- 
mate that the conversion from Sloan AB mags to flux density 
units introduces an uncertainty in the colour fgi / fii of ~10 
percent. It should also be noted that the residual superhump 
will introduce an additional uncertainty in the colour of the 
flare. Since the superhump originates in the outer parts of 
the disc, it is most likely to have a cool spectrum. Theo- 
retical models for the behaviour of the quiescent accretion 
disc in V404 Cyg indicates that the outer disc is optically 
thick and has a temperature of ~2000 K (Narayan, Barret 
& McClintock 1997). If this also represents the colour of 
the superhump, then assuming the superhump radiates as a 
blackbody, we estimate that it would add an additional un- 
certainty of 2 percent in the determination of fgi /fii . Note 
that these uncertainties are much less than that due to the 
reddening (see section l^TTl . 



4.1 The characteristics of the flares 

Although large amplitude flares are observed on more than 
one night, it is only on Sep 10 that a large flare is observed 
throughout the event. The shape of the large flares seem 
to be symmetric with similar rise and decay times. Super- 
imposed on these large flares are many rapid events which 
last typically ~0.5 hr. Furthermore, superimposed on these 
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Figure 3. Detailed plots of some individual distinct fiares in the 
i' -band lightcurve. The labelling refers to the events in Figure 2 
and Table 2. 



flares, are still shorter term events on timescales of mins 
which show unresolved peaks (see Figure 3). This makes it 
difficult to determine the exact shape of the fiare events; 
there is an indication that some of the fiares on the short- 
est timescales have a burst-like profile [e.g. see flare (b) on 
Sep 13]. However, in general all the flare events appear to 
have a linear rise and decay. In Table|5|we estimate some of 
the characteristics of the individual distinct flares which are 
shown in Figure 3. The flares typically have a rise time of 
2 mins and a duration of 6 mins. 

At first sight, once clearly notices that the monochro- 
matic flux in i' is larger than that it g' , at the peak of the 
large flares, i.e the flare at its peak has a relatively red spec- 
trum. In general the colour flux ratio fg' / fi^ of these flares 
increases, becoming more blue at the peak. After the peak is 
reached, the flux and colour decreases becoming more red. 
The exception is the large flare on Sep 13; the colour seems to 
increase continuously as the flare proceeds. The most promi- 
nent and distinct rapid flares are marked in Figure 2. The 
mean flux ratio fgi /fit of the large flares is ~0.4. If we use 
the extreme values for the reddening ^v=4.0 and Ay=2.2, 
then the fiux ratio fgi j fii changes by a factor of 0.53 and 
1.36 respectively. 

To test if the flaring in the g' and i! bands are simul- 
taneous we calculated the cross-correlation function. The 
cross-correlation function was symmetric and no time lag 
was found, so we conclude that the flaring in the two bands 
is simultaneous (see Figure 4). 



5 THE POWER DENSITY SPECTRUM 

To compute the PDS of the "slow" data, we use the same 
flare lightcurves as described in section^] but since we are 
only interesting in sub-orbital variability, and in order to 
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Figure 4. Cross-correlation function between the g' and i' -band 
lightcurves. A postive lag would indicate that the g' variations lag 
behind those in i' . 



Table 2. Properties of the V404Cyg flare events. The rise and 
decay times are obtained from linear flts to the i' band flare pro- 
file. 



Date 


Flare 


Rise 


Decay 


PDS power-law index* 




No. 


(s) 


(s) 


g' i' 


Sep 10 


al 


240 


110 


-1.23 -1.56 




bi 


104 


100 






b2 


57 


58 






ci 


65 


55 




Sep 12 


al 


322 


242 


-0.96 -1.32 




bi 


62 


54 






b2 


88 


44 






b3 


62 


187 






ci 


175 


241 




Sep 13 


al 


170 


136 


-1.17 -1.49 




bi 


20 


96 






b2 


22 


78 






ci 


90 


130 






c2 


140 


140 




Power-law index 1- 


(T errors 


of 0.03. 





preserve the flux in each band, we detrend the data for each 
night using the nightly mean and add on the mean flux level 
of all the data. Although the samphng of ULTRACAMis 
perfectly uniform, we use the method of Lomb-Scargle to 
compute the periodograms (Press et al. 1992). The Lomb- 
Scargle method was chosen so that we could combine consec- 
utive runs and confidently compute the periodogram. Also, 
in order to allow a direct comparison with X-ray PDS, we 
use the same normalisation method as is commonly used 
in X-ray astronomy, where the power is normalised to the 
fractional root mean amplitude squared per hertz. The ad- 
vantage of this kind of normalisation is that fractional rms 
ampli tudes can be direc tly estimated from the level of the 
PDS (Ivan der Klislll995l) . 



To compute the PDS we use the constraints imposed 
by the Nyquist frequency and the typical duration of each 
observation. We then average the PDS for each observing 
sequence. We bin and fit th e PDS in logarithmic space 
jPaoadakis fc Lawrencelll993h and the errors in each bin are 
determined from the standard deviation of the points within 
each bin. The white noise level was subtracted by fitting the 
highest frequencies with a white noise plus red noise model. 



5.1 The PDS of the "slow" data 

To determine the PDS for the "slow" data we only use the 
data with the longest baseline i.e. during the nights of Sep 
10 and 12. This allows us to extend the PDS to the lowest 
possible frequencies. The derived PDS is shown in Figure 5. 
The shape of the g' and i' band PDS are well described by 
a power-law with slopes of -1.24±0.02 and -1.37±0.02 re- 
spectively (l-(7 errors). On each individual night, the slope 
of the power- law changes (see Tablel^J. To examine the sys- 
tematic effects in our analysis, we also computed the PDS 
of the comparison star. The PDS was found to be flat, ex- 
actly as expected for white noise data. Furthermore, the 
auto-correlation function of the comparison star is a delta 
function, implying that the data are not correlated. A QPO 
feature with a centroid frequency at 0.78 mHz {= 21.5 mins) 
is present, not surprising given that one can clearly see vari- 
ability on timescales of 0.5-1.0 hrs in the lightcurves. The 
fractional root mean squared amplitude is 3 and 1.3 percent 
for the g' and i' band QPOs and was computed by integrat- 
ing the Gaussian function and taking the square root. 



5.2 The PDS of the "fast" data 

We used runl4(a), runl4(d) and runl4(e) to compute the 
PDS. The other runs were not included because of the high 
noise level of the data, the result of an increase in the elec- 
tronic bias level (see section|5J. The computed "fast" PDS 
is shown in Figure 5. The general shape of the PDS has a 
power-law shape with a slope of — 1.16±0.04 (l-cr errors). 
As a check we also computed the PDS of the comparison 
star. Again the derived PDS was flat, which is exactly as 
expected for white noise data (see Figure 5). 



6 DISCUSSION 

6.1 The slope of the PDS 

The outer regions of the accretion discs in quiescent SXTs 
and dwarf novae (DNe) are assumed to be very similar. 
Therefore if the optical variability observed in the quiescent 
SXTs originates from these outer regions, then they should 
have si milar properties to those observed in quiescent DNe. 
lBruc3 |^w3) has performed a study of the rapid variabil- 
ity in cataclysmic variables. The slope of the PDS in quies- 
cent DNe (-1.6 to -2.6) seems to be in general steeper than 
those observed in quiescent SXTs [-1.0 to -1.52: IZurita et aH 
(2003) and Hyncs ct al. (2003) ]. This suggests that there is 
either more low-frequency variability in the lightcurves of 
quiescent DNe compared to quiescent SXTs, or that there is 
more high-frequency variability in the SXTs. Note that the 
former could be attributed to intrinsic variability associated 
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Figure 5. The power density spectrum of V404Cyg. Prom top to 
bottom: V404 Cyg 3' and i' "slow" PDS, V404Cyg i' "fast" PDS 
and tlie comparison star PDS. 

with the stream/disc impact region or the white dwarf (e.g. 
IWarnej ri995). 

On average the slope of the i' PDS in V404 Cyg is 
steeper than the g' slope, suggesting that there is more low- 
frequency disc variability in i' compared to g' . This differ- 
ence could be attributed to more variability from outer re- 
gions of the disc, such as the stream/disc impact region or 
more high frequency variability in g' from the inner disk. 
Thus one expects the slope of the PDS to be dependent on 
the observed waveband . 

iHvnes et al.l J2003D determined the PDS of the quies- 
cent SXT A0620-00 and found a possible break in the PDS 
at 0.95 mHz, a fe ature which is also seen in the low/hard 
state of SXTs fsee lWiinands fc van der Kli!j|l999l and refer- 
ences therein) . Although the origin of the break frequency is 
uncertain, it could be related to the size of the inner advec- 
tive region; a low break frequency in quiescence would then 
be expected, as the advective region should be larger. If the 
break frequency is associated with a dynamical timescale, 
then it should vary as R}'^ ^ where R is the radius of the 
ad vective region (see section lF!^ . Furthermore, as noted 
bv iHvnes et alj (|200^ the break frequency most probably 
scales with some characteristic length scale, which could de- 
pend on other parameters, as well as the inner disc radius. 

6.2 The QPO feature and the ADAF model 

iNaravan. McClintock fc lll996ll and iNaravan et alJ 
il997il showed that the observations of quiescent BH 
SXTs can be explained by a two-component accretion 
flow model. The geometry of the flow consists of an inner 
hot advection-dominated accretion flow (ADAF) that 
extends from the black hole horizon to a transition radius 
rtr=10"'~^ rsch {rsch.— GMjf? is the Schwarzschild radius 
for a black hole with mass M) and a thin accretion disk 
that extends from rtr to an outer radius rout i.e. the radius 
of the accretion disc. A key feature of the ADAF model of 



quiescent BH SXTs is its low radiative efflciency. The bulk 
of the viscously dissipated energy is stored in the gas and is 
advected with the flow into the black hole. This naturally 
explains the unusually low luminosity of BH SXTs, since 
in NS SXTs, all the advected energy is radiated from the 
neutron star surfa ce, resulting in a m uch higher overall 
radiative efficiency (INaravan fc Yi|[i995^ . 

An ADAF has turbulent gas at all radii, with a variety 
of timescales, ranging from a slow timescale at the transi- 
tion radius down to nearly the free-fall time close to the 
black hole. In principle, interactions between the hot inner 
ADAF and the cool, outer thin disk, at or near the tran- 
sition radius, can be a source of optical variability, due to 
synchrotron emission by the hot electrons in the ADAF. 
For an ADAF the variability could be quasi-periodic and 
would have a characteristic timescale given by a multiple 
of the Keplerian rotation period at rtr tK ~ 2ttR/vk = 
27r(GM/i?^)-^''2 ~ (M/10Mo)(r/100)^ '^s. where R is the 
absolute radius and r is in units of rsch. One also expects 
slower variat ions, since the mass supply to the ADAF orig- 
inates at rt,. iNaravan et al]|l997^ . It should be noted that 
it is difficult to produce such rapi d variations using the thin 
disc models. INaravan et alJ^1997^ used the maximum veloc- 
ity observed in the Ha emission line arising from the quies- 
cent thin accretion disc to estimate rtr- However, they note 
that measuring the maximum velocity is difficult, and so 
their estimate of 10'* '*rsch for rtr should only be considered 
as an upper limit. If we assume that the 0.78 mHz QPO fea- 
ture observed in the PDS (see section|^ is the dynamical 
timescale at the transition radius, then rtr li es at 10^ °rsch. 
which is consistent with that estimated by INaravan et alJ 
(1997J. We can write rtr in terms of the accretion disc ra- 
dius. For the binary parameters of V404 Cyg the outer radius 
of the disc lies at rout=2.3 x lO^rsch or 12 Rq, and so rtr 
lies at 0.04 7?discOr 0.5 Rq. 

6.3 Comparison of the V404 Cyg QPO features 
with 4U 1626-67 

It is interesting to note the similar timescales of the QPOs 
we have discovered in V404 Cyg with those found in vari- 
ous X-ray /UV/optica l observations of the double-degenerate 
LMXB, 4U 1626-67. IChakrabartv et alj J200ll) used fast 
CCD photometry to reveal the presence of 3 oscillations: the 
130 mHz pulsation (i.e. the well-known 7.66 s X-ray/optical 
pulsation of the neutron star), a 48 mHz (~20s) QPO (a 
broad oscillation seen before in X-ray, optical and UV wave- 
bands) and a ImHz QPO. This latter oscillation had not 
been seen before, is not present in any X-ray observations, is 
quasi-sinusoidal in ap pearance, and is dramatica lly stronger 
in far UV HST data. IChakrabartv et al.l J200ll) concluded 
that the first 2 oscillations were clearly X-ray driven and 
the optical/UV analogues were due to reprocessing in the 
disc/companion star. 

However, the 1 mHz QPO origin is extremely puzzling. 
The timescale is actually the same as the pronounced X-ray 
and optical flaring seen from 4U1626-67 during the 1970s 
and 80s, but this flaring behaviour has vanished since the 
torque reversal that occurre d in 1990. With no X-ray oscil- 
lation to drive this process, IChakrabartv et "al] i200 j) sug- 
gested that the 1 mHz QPO is due to "warping" of the in- 
ner accretion disc and is of an appropriate timescale for the 
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radius of the inner disc that is set by the magnetospheric 
corotation radius. In this way, the modulation would be ex- 
pected to be much greater in the far UV (inner disc) than 
the optical (in the outer, shielded regions). It is possible 
that, in spite of the fact that V404Cyg almost certainly 
hosts a ~12 M0 black hole, this analogy can be applied 
from 4U 1626-67. The difference is in the origin of the trun- 
cated inner disc, which for V404 Cy g in X-ray q uiescence is 
expected to be in an ADAF phase ijNaravan et aljilQQTi '). 

6.4 The emission mechanism of the flares 

As one can see from Figure 2, on most of the nights there are 
numerous rapid flares (<0.5hr) superimposed on the single 
large event (~few hrs). Given the uncertainties in the data, 
it is difficult to determine the shape and hence the physical 
parameters of these rapid events. We can however comment 
on the general physical properties of the flares using the large 
events (~few hrs). Using the flare flux density lightcurves we 
determine the flux density ratio fgi / fii and the blackbody 
temperature Tbb and equivalent radius Rbb- To determine 
the colour temperature corresponding to a given flux den- 
sity ratio, we integrate blackbody functions with the CCD 
and Sloan filter response functions and then determine the 
fg' / fi' flux density ratio. Given this blackbody temperature 
we can then determine the corresponding radius of the region 
that pr oduces the observed dereddened flux at a distance of 

3.5 kpc dCasares et alJl993^ . Also, using the flux density ra- 
tio, the power-law index a of the spectrum {Fv oi v") can 
be obtained. One can easily show that the power-law index 
is given by a = A log f./A log j/ = 4.83 log(/<,,//iO. 

For Sep 9, the flux density ratio /g/ / f^i remains constant 
at 0.5; a and Tbb have constant values of ~-2 and 5100 K 
respectively. For the single large flare events on Sep 10 and 
13, fgi/fi' is correlated with the flux of the flare; fg'/fi' 
peaks at ~0.3; a and Tbb peaks at ~-2.5 and 4100 K re- 
spectively. On Sep 12 one can see that the blackbody radius 
and flare flux are correlated i.e. as the flare flux increases 
and decreases so does the equivalent blackbody radius. The 
blackbody radius increases from ~lR0to ~2R0. The same 
is seen during the decay of the large flare on Sep 10, but 
to a lesser extent. However, it should be noted that the un- 
certainties in the reddening (see section|^, which typically 
correspond to ±700 K in the Tbb estimates, are such that 
we cannot comment on any colour or temperature variation 
during the flares. 

The peak flux produced by the large flare arises from a 
region that has a projected blackbody radius of 2 Rq and 
the rapid flares most likely from smaller regions. In other 
words, the emitting area at the peak of the large flare covers 
at least 3 percent of the total disc emitting area. Note that 
given the data presented here, we cannot state where in the 
disc these large flares arise. However, it is most likely that 
the large flares originate in the outer regions of the disc (see 
section l^3l . 

During the large flares observed in V404 Cyg the power- 
law index of the spectral energy distribution is ~-2.0 (a 
range of -1.5 to -2.5, depending on the flare). It is difficult 
to explain this in terms of optically thin synchrotron emis- 
sion, unless the electrons follow a much steeper power-law 
energy distribution {N{E)dE oc E^dE). For solar and stellar 
flares the power-law index of the electron energy distribu- 



tion p is ~-2 iCrosbv. Aschwanden fc Schmitt|ll993l) . which 
corresponds to a frequency spectrum with a power-law index 
of q:=-0.5 {Fi, oc r/'-^^'''/^). The frequency power-law index 
oberved in V404 Cyg a of ~-2.0 implies a much steeper in- 
dex for the electron energy distribution of p=-5. However, 
given the uncertainites in the reddening and therefore in the 
flux density ratio fgi / fii (see section|3J , we find a range of 
values for p of -3.6 to -7.5. It could be that the electrons 
are thermal and thus do not follow a power-law energy dis- 
tribution. Also, given that the geometry and physics around 
a black hole is presumably more extreme, such a steep in- 
dex for the electron energy distribution is not completely 
implausible. 

The simple model for the flares of blackbody radiation 
from a heated region of the disc's photosphere (i.e. photo- 
spheric in origin) is unlikely. A more likely model for a ther- 
mal flare is emission from an optically thin layer of recom- 
bining hydrogen, which is essentially the mechanism gen- 
erally accepted for solar flares. We therefore determine the 
continuum emission spectrum of an LTE slab of hydrogen 
and then calculate the expected flux density ratio fgi / fii in 
the g' and i' bands using the same method as the blackbody 
case. We find that a continuum slab temperature of 8000 K 
with a baryon column de nsity of lO^'^ cm~^ (taken from 
IWelsh. Home fc Okelll99l) gives /g'//»'=0.30; propagating 
the uncertainty in the reddening (see section|3J, we find the 
range for the slab temperature to be 5500 K and 9100 K. It 
should be noted that the hydrogen slab model does not in- 
clude the contribution from Balmer emission lines. Since the 
g' filter includes H/3, our determination of the flux density 
ratio fgi I fii is biased. However, we can estimate the amount 
of bias given the H/3 flux and equivale nt width, 0.20 mjy and 
lOA respectively dCasares et alJll993ll . Since the g' filter has 
a width of ~2000A the contribution of the H/3 emission line 
to the g flux density in V404 Cyg is only 1 percent, in con- 
strast to that in AE Aqr JWelsh. Home fc Okelll993l) . 

Since we only have observations in two wavebands, we 
cannot accurately determine the physical parameters of the 
flares. However, it should be noted that the spectrum of 
the flare seems to be simil ar to that observed in AEAqr 
iWelsh. Home fc OkdlT99l . Using the colour information 
for the large flare events, we can say that the large flares 
arise from a region which extends a projected area which 
has an equivalent i?BB of ~2R0. One should regard the 
equivalent -Rbb estimate as a lower limit, because the emis- 
sion mechanism is unlikely to be blackbody. In general the 
rapid flares seem to arise from optically thin gas which has 
a peak temperature of ~8000 K in the continuum. However, 
to obtain a more accurate estimate of the physical condi- 
tions of the flares, one must resolve the Balmer jump and 
Paschen continuum. 

Spectroscopic observations suggest that the whole disc 
participates in the large flar es, since the whole Hq emis- 
sion line profile changes llH^nes et alJl2002l) . This seems to 
contradict our earlier results, which suggest that the large 
flares arise from regions that cover only 3 percent of the 
disc's total area. However, this most likely indicates that 
the emission is optically thin, or that the emission area is 
small. Also, given the low resolution of the spectroscopy, we 
can only really say that material from all around the disc is 
contributing, given that the line proflle changes observed ex- 
tend across the whole line profile and appear double-peaked. 
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It is not possible to exclude the flares from arising from only 
the inner regions of the disc. 

6.5 Where are the flares produced? 

IZurita et al.l i2003l) have presented arguments that suggest 
the most likely mechanism for the flares is by local mag- 
netic reconnection events in or above the accretion disc. It 
is believed that a dynamo mechanism, driven by the strong 
shear prod uced by differential r otation, operates in accre- 
tion discs jTout fc Pringldll992^ . Regions of oppositely di- 
rected magnetic fields develop within the disc or between 
the d isc and corona and reconnect exp losively, resulting in a 
flare jHaswelL Taiima fc Sakaill992bll . The duration of the 
flare is di rectly related to th e shearing timescale. As sug- 
gested bv IZurita et al.l i2003tl . assuming that the range of 
flare durations observed in the SXTs reflect their location 
in the accretion disc, then the flares with the longest dura- 
tion represent the shearing timescale in the outer parts of the 
disc. Thus the maximum flare duration should be limited by 
the outer disc radius and V404 Cyg should show short-term 
variability, similar to what is observed in A0620-00. The 
very rapid flares (~mins) observed in V404Cyg (see Fig- 
ure 3) have a similar timescale to those observed in A0620- 
00 iHvncs ct al. 2003) and seems to support the scenario 
that the large flares (~ few hrs) are produced in regions fur- 
ther out in the disc than the rapid (~0.5 hr) and very rapid 
(~mins) flares. 

Using the colour information about the large events 
(~hrs), we can say that the large flares are produced from 
a region which extends an equivalent blackbody radius of 
~2R0. We cannot determine the emitting area of the rapid 
flares (~0.5hr) as the result may be biased due to the large 
flares. However, if our interpretation of the 21.5 min QPO 
feature is correct, then these rapid flares are most likely pro- 
duced near the transition radius of the ADAF, in the inner 
parts of the disc. We cannot determine the exact mechanism 
for the flare production. The flares could arise in optically 
thin gas with a temperature of ~8000 K, most likely in a 
corona above the accretion disc, or from optically thin syn- 
chr otron emission. 

iNavakshin &: SvenssonI ||20^) have argued that an 
ADAF model is not required to reproduce the observed X- 
ray luminosities of quiescent SXTs. They propose a model 
which consists of a standard thin disc with a hot corona, 
powered by magnetic flares. Although this scenario seems 
appealing, it is difficult to see how the thin disc model 
could produce rapid variations. However, it should be noted 
that the explosive magnetic recon nection model described in 
iHaswell. Taiima fc Sakail lll992bft could produce rapid varia- 
tions from a thin disc, similar to an ADAF. Our observations 
suggest that for the rapid flares ('^O.Shr) at least, they are 
produced near the transition radius. However, it is only by 
obtaining precise velocity information about the flares, can 
we begin to place constraints on the exact location of the 
flares in the accretion disc. 
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